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Abstract 

In this paper we use a generalization of Oevel's theorem about master symmetries 
to relate them with superintegrability and quadratic algebras. 



I Introduction 

In this article a general framework is built in terms of master symmetries and recursion 
operator to provide superintegrability and quadratic algebras. 

This is applied to the isotropic harmonic oscillator, the rational Calogero-Moser system 
and the "Goldfish" model. 



II Background 

Let M be a differentiable (C°°) manifold of finite dimension and A a bivector (a 2-times 
contravariant skew-symmetric tensor field) on M. Associated with A there is a natural 
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morphism A* from the cotangent bundle T*M into the tangent bundle TM defined, for all 

a,(3e T*M, by 

< A"(a),/5 >= A(a,/5), (2.1) 

where < ., . > denotes the usual coupling between 1-forms and vector fields. 
We also define a bilinear map from C°°(M) x C°°(M) into C°°(M) by 

{f,g} = A{df,dg), f,geC^{M). (2.2) 
Due to the properties of A, this bracket satisfies 

PBl {f,g} = -{g,f} skew- symmetry 

PB2 {fh, g} = f{h, g} + {/, h}g Leibniz rule 

We say that (M, A) is a Poisson manifold, and A is a Poisson tensor, if, in addition, the 
bracket ( p.2|) satisfies the Jacobi identity 

PBS {/, {g, h}} + {g, {h, /}} + {h, {/, g}} = 0. 

which is equivalent to the vanishing of the Schouten-Nijenhuis bracket [A, A]. 
We call a vector field X an infinitesimal Poisson automorphism if it satisfies 

X{f,g} = {X{f),g} + {f,X{g)}, f,g G C^{M). 

Given a differentiable function H on M the Hamiltonian vector field associated with H 
is the vector field defined by 

Xh{x) = A^{x){dH{x)), xeM. (2.3) 

If A is a Poisson tensor then 

[Xf,X,] = [A^{df),A^{dg)] = A\d{f,g}) = X^j^.y (2.4) 

that is. A" is a Lie algebra homomorphism between the Lie algebra of differentiable functions 
(C°°(M), {.,.}) and the Lie algebra of vector fields {A{M), [.,.]). 

An integral of motion of H, or of X^, is a differentiable function F such that 



{H,F} = Xh{F)=0, 
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(2.5) 



so is constant along the orbits of the Hamiltonian system Xh 

f = ^h(0). (2.6) 

By a Nijenhuis operator i? in a manifold M we mean a (1, 1)— tensor satifying, for all 
vector fields Z in M, 

Cr^z)R = RCzR. (2.7) 

The Nijenhuis operators transform closed 1-forms into closed 1-forms in the following 
sense 

Proposition 1 Let R be a Nijenhuis operator and a a closed 1—form such that ai = *Ra 
is also closed then, for all i G N, = (where means i*^ iterates of^R) are closed. 

By a conformal vector field of a tensor W we mean a vector field Z such that CzW = cW, 
for some constant c G IR. 



With Oevel's Theorem |11|, a Nijenhuis operators helps to define new symmetries if a 



conformal vector field is known. 

Theorem 2 (Oevel) Let R be a recursion operator of Xq G A{M), this means CxR = 0, 
and Zq G A{M) a conformal vector field of Xq and R such that 

^Zo^o = -^-^05 ^ZqR = fJ'R-, A, /i G H. 

If R is also a Nijenhuis operator then, defining the sequences Xn = R^Xq and Zn = R^Zq, 
n G IN, itie have, for all n,m & No 

Cz,.R = fiR^+\ 

and 

[Zn, XJ^ = (A + m^)Xn+m- 

The Zj's are called master symmetries or symmetries of second order oi the vector fields 
Xj because [[Zj, X^], X^] = but [Zj, Xj\ ^ and they help us define new symmetries of the 
system. We call master integrals to functions G which may not be integrals of motion of the 
system but they induce an integral X{G). 
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Now let us consider two linearly independent Poisson structures Aq and Ai in M. We say 
that they are compatible if Aq + Ai is again a Poisson tensor. The compatibility condition 
is equivalent to the vanishing of the Schouten-Nijenhuis bracket [Aq, Ai]. 

A vector field X G A{M) is said to be hihamiltonian if it is Hamiltonian with respect 
to two independent compatible Poisson tensors, that is, if there exist two functions H, F & 
C^{M) such that 

X = Kl{dH) = K{{dF). (2.8) 

There is an important example of bihamiltonian manifold: the Poisson-Nijenhuis man- 
ifold [^]. This manifold is special because one of the Poisson tensors is obtain from the 
other by means of a Nijenhuis operator. For instance, a bihamiltonian manifold (M, Aq, Ai) 
such that the Poisson tensor Aq is non-singular (symplectic manifold), is a Poisson-Nijenhuis 
manifold with Nijenhuis operator R = AJAq"^. If Xi = Ao((iifi) = Af (rfi/o) is a bihamilto- 
nian system then we may define a sequence of symmetries of X, Xi = R^^^Xi and, if the 
first cohomology group is trivial, a sequence of integrals of motion in involution (//j)jgN such 
that X, = kl{dHi), ie dHi = ^R\dHo), i G IN. 

Ill The superintegrability and the cubic algebra 

Let X be a vector field on a manifold M of dimension n. It is called maximally superintegrable 
if it possesses n — 1 functionally independent first integrals. 

There are several examples of maximally superintegrable systems, some of them are 
shown in the last section. 

The superintegrability may be a consequence of the existence of a sufficient number of 
master integrals, as the next proposition shows. 

Proposition 3 Let X be a vector field on a manifold M and G.F & C°°{M) master integrals 
of X . Then the function 

L = X{G)F - X{F)G (3.9) 

is an integral of the vector field X . 
Proof: 

Just notice that 

X{L) = X{X{G))F + X{G)X{F)-X{X{F))G-X{F)X{G) 
= X{X{G))F-X{X{F))G = 0, 
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because that X{F) and X{G) are integrals of the system. 



Remcirk 4 If X is a vector field on a manifold of dimension 2n and if n functionally in- 
dependent master integrals are known Gi, . . . ,Gn then we can define the integrals of motion 
Fi = X{Gi) and Lij = X{Gi)Gj — GiX{Gj) which may provide the superintegrahility of the 
system if 2n — 1 of them are functionally independent. 

Theorem 5 Let X be a Poisson infinitesimal automorphism on a Poisson manifold (M, {., .}) 
Suppose there exist master integrals of X, Gi, such that {Gi,X{Gi), i & J C. INT} is a basis 
of a Lie subalgebra of (C°°(M), {.,.}) 

Then, for all i, j e N the functions X{Gi) and LGi,Gj = X{Gi)Gj — X{Gj)Gi generate a 
cubic algebra for the Poison bracket. 

Proof: 

Once {Gi,X{Gi),i e J} generates a Lie algebra, for each i,j,k e J there exist constants 
a,^,-, b^ A such that 

{G,, G,-} = El^t-G/c + hl^X{G^)). (3.10) 
fceJ 

Applying X twice to the last equation we obtain 

X({X(GO, (G,-)} + {G,,X(G,)}) = at.X(G,)), 

fceJ 

so 

{X(G,),X(G,)} = 0. 

Writing {Gj, X{Gf)) = ^{c'^jX{Gk) + d'^jG^), with c^j, d^j constants, and noticing that 

keJ 

X{{G,,X{Gj)}) = {X{G,),X{Gj)} = 0, 

we have ^ d\^^X{Gk) = 0, which yields d^^ = 0. 
fceJ 

Thus 

{LG„G,,X{Gk)} = X(G,){G,-,X(Gfc)}-X(G,){G,,X(Gfe)} 
= Y.KkX{GmGi) - 4X(G,)X(G,)] 
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and 

{LGi,Gj,LGk,Gh} = 

= X{Gk)[L{Gi,Gh},G, - L{Gj,Gh},Gi] + X{Gh)[L{Gj,Gk},Gi " L{Gi,Gk},Gj] 

+LG^,GdG^,XiG,)} + LG,,GAG„XiG,)} 
-LG^,Gh{Gi,X{Gk)} - LGi,Gk{Gj,X{Gh)} 
= Yl\.'^iGk){a[hLGi,Gj + b\ f^Lx(Gi),Gj - a!'j f^LGi,Gi - h] f^Lx(Gi),Gi) 

+X{Gh){a}j,kLGuGr + b^j,kLx(G,)r,, - (^[.kLGiXi, - b\}.Lx{G,).G,) 
+4,hLG„G,X{Gi) + c)^kLG.,G,X{Gi) - c{kLG„G,X{Gi) - 4,Lg„g,X(G,)] 

= Yl^iGk){a\hLGi,G:j - a] f^LGi,Gi) + X{Gh){a],.LGi,Gi - a\,,LGi,Gj) 

I 

+cl,hLGj,GkX{Gi) + c] j^LGi,GhX{Gi) - c\ f,LG^,GhX{Gi) - c] f^LGi,GkX{Gi) 
Hh\,HX{G,) - h]^^X{G,))X{G{)X{Gk) + {h]^kX{G,) - h{,X{G,))X{G{)X{Gn)\ 

This last expression being cubic in the quantities X{Gi) and Lg^^Gji we refer to this result 
as that they generate a cubic algebra. 

■ 

Corollary 6 If for each i,j E N, J2ibijX{Gi) = hiX{Gj) — bjX{Gi) or all the constants 
b^j are zero then the integrals of motion genereate a quadratic algebra. 

Remctrk 7 Notice that in the above proposition we could have just demanded that d{Gi, Gj} 
be a linear combination of the dG 's and the dX{G) 's. 

Theorem 8 (Generalization of Oevel's theorem) LetX be a vector field on a manifold 
M, R a Nijenhuis operator which is also a recursion operator of X and P a (1, 1) — tensor 
satisfying 

CxP = a{R) (3.11) 

and 

jCpxR^b(R), (3.12) 

with a{R), b{R) polynomials in R. Then, defining the sequences Xi = R^X, Yi = R\PX), 
i e INo, we have 

[X„X,] = 0, (3.13) 
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[X„ Y,] = a{R){X,+,) - ih{R){X,^,.i) (3.14) 

[F„y,] = (j-z)f>(i?)l"i+,-i. (3.15) 

Proof: 

The proof is similar to the original Oevel's theorem proof, which can be seen at ||lT| . ■ 



Suppose that M has trivial first cohomology group and is endowed with a non-degene- 
rated Poisson structure A such that -RA" = A"*/? (this means that the tensor RK is a bivec- 
tor). Futhermore suppose that there exist functions such that X = A^[dHi) = i?A^((iifo) 
and Y = A\dGi) = RA^{dGo). 

Then Proposition |l| ensures us that the 1— forms 

ai = 'R\dHi), Pi = 'R'idGi), i G IN 

are closed and we can consider them exact because of the triviality of the first cohomology 
group. 

Write ai = dHi and I3i = dOi, for alH G IN. 
First notice that RA being a bivector yields 

X,{Hj) = <X,,dHj>=<X,+„dH^> 
= R'+^A\dHi,dHi) = 0. 



Moreover (|3.14|) ensures that the G's are master integrals of the X's because 



X,(X,(G,)) = X,i{Hi,Gj}) = -[X,,Y,]iHi) 

= {ib{R)X,+j_^-a{R)Xi+^){H,) = Q, 



relation ( |3.15| ) implies 

d{G,,G^} = {j - i)h{'R)dGi+,_^ 



and 



{X,(G',),G,} = -y,(X,(G,)) = [X„n](G,) - X,{Yk{G^)) 

= ib{R)X,+k-i{Gj) - a{R)Xi+k{G,) - X,({Gfc, G,}) 

= ib{R)X,{G,+k-i) - aiR)Xi{Gj+k) - (j - k)b{R)Xi{Gk+,-i) 

= {z + k-j)b{R)X,{G,+,-i)-a{R)X,{G,+k). 



So {Xi{Gj),Gk} can be written as a linear combination of the Xi{G)'s. 

Now we can apply Theorem ^ with the 6's equal to zero, and guarantee that, for each 
2 G IN the integrals of Xj, Xi{Gj) and L^j = Xi{Gk)Gj — Xi{Gj)Gk, j, k G INq, close quadrat- 
ically under the Poisson bracket. 

Futhermore notice that 

[RA, A](A«-ir) = Cry A + (>CyA) o'R+ [CyR) o A 

so, as Y and RY are Hamiltonian vector fields, we have 

[RA,A]{A^-^Y) = (CyR) oA = b{R)A. 

Thus, if RA is a Poisson tensor then it is compatible with A, because i? is a Nijenhuis 
operator, and b{R)A = 0. 
But this implies that 

b{R)Xi = b{R)A\dHi) = and b{R)Yi = b{R)A\dGi) = 0, 

so the relations ( |3.13| ), ( p.l4| ) and ( p.l5| ) become 

[X„ X,] = [y„ Y,] = 0; [Xi, Y^] = a{R)X,^, 

and 

d{G,,Gi] = 0; d{Gu,X,{G,)} = a{R)X,{G,+k). 



Although in the last procedure we need two master integrals, Hamiltonians of the master 
symmetries, to generate a new sequence of integrals of motion, we may construct it only 
knowing one master integral of all the vector fields. 

Proposition 9 Under the conditions of proposition^, suppose there exists a master integral 
G of all the vector fields Xi, i G Nq, then the functions Gi = Yi{G), i G INq are also master 
integrals of the same vector fields and, for each k G INq, 

Ll^ = Xk{G,)Gj - Xk{Gj)G„ for all i,j G INq, 

are integrals of Xj.. 
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Proof: 

Due to relation (|3.14|) , we have, for all i,k E INq 



Xk{G,) = [Xk,Yi\iG) + Y,iXk{G)) 

= a{R)X,+k{G) - kb{K)X,+k-iiG) + Y,{Xk{G)). 

But 

X,{Y,{Xk{G))) = a{R)X,+k{Xk{G)) - kh{R)X,+k-i{Xk{G)) + y,(X,(Xfe(G'))) = 0, 

because Xk{G) is an integral of all the Xj, so Yi{Xk{G)) is an integral of Xk and Xk{Gi) is 
then an integral of X^. Thus Gi is a master integral of all the X^. 



Remark 10 It seems that in the previous procedure only one Hamiltonian Gq would be 
necessary, but note that the procedure applied to Gq yields Gi = Yi{Go) = R^A{dGo, dGo) = 0, 
i.e. all master integrals are zero. 

IV Examples 

Example 11 (Isotropic Harmonic Oscillator and the Fernandes' Theorem) 

The Isotropic Harmonic Oscillator is the Hamiltonian system in (JR^"" \ (qi,pi)), defined 
by the Hamiltonian 

n 1 

^ = E^b? + a- (4-16) 

i=i ^ 

It is well known that this system is completely integrable with constants of motion in 
involution 

Ei = ^ip^ + q^), ^ = l,...,n, (4.17) 
and it is bi-Hamiltonian with respect to the compatible Poisson tensors 

^ A A — 
:i dpi dq, 



" (9 8 



and 
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The Hamiltonian vector field can be expressed as 

= t.iP^^ - = ^Udn) = AlidHo), (4.20) 

i=l ^Qi 

with Hq = ln(i?i) + . . . + ln(£'„). Defining the recursion operator of the system 

R = A\aI' =J2Ei{—(g)dq, + — (g) dp^), (4.21) 

i=l OPi 

a sequence of Hamiltonians and of Hamihonian vector fields can be defined 

Xi = N'-^Xn = kl{dH,), i = l,...,n. (4.22) 
Notice that the (1, l)-tensor 

P = ll ® dqi + 7^- ® dpi), (4.23) 

j=l (jpi 

with Lpi = arcsin( ), satisfies the conditions of Theorem § with a{R) = Id and b{R) = 

0. So defining the sequence 



we have 



and also 



n = R'PX = Y: V^EHP^T- - q^T-), (4-24) 
i—l OQi Opi 



[Yi,Xj] = -Xi+j 



dZ,{Hj) = 



[y„i?^Ao] = -i?*+^Ao 

i.e., the vector fields Yi are not Hamiltonian with respect to any of the Poisson structures. 

n 

But let us define the functions G = ^ E^ipk and notice that 

2=1 

X,{G) = Y: El-\pk^ - qu^){G) = iH,. (4.25) 
^ oqk opk 

So the functions Gj = Yi{G) = J2k VkE]^'^ satisfy 

X,{G,)=X;^,{G) = {i+3)H,+,, 
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= {i + j)Xi^,{Gk) = {i + ])X,+j+k{G) 
= ii + j + k)Hi+j+k 

and 

{Gi,Gj} = Ao{dG,,dGj) = Ao{dZ,{G),dZ,{G)) 

n n 

= Ao(rf(E i^^+Vfc), ^(E ^rV/.)) = - J)Z^+,iG) = (j - 

k=l h=l 

Therefore, for each k G Nq the integrals of motions of Xk, Xk{Gj)'s, and L^j = 
Xk{Gi)Gj — Xk{Gj)Gi = {i + k)Hi^kGj — (j + k)Hj^kGi dose quadratically under the 
Poisson bracket defined by Aq. 

Now let us consider a little more general configuration, in which the isotropic harmonic 
oscillator is a particular case. 

Given a completely integrable Hamiltonian system (M^", u;, H) in a symplectic manifold, 
Fernandes establishes necessary and sufficient conditions for the existence of a second 
Poisson structure giving the complete integrability of the system, in a neighborhood of a 
fixed invariant torus. Without lost of generality let us consider 

(M^" = 1R" X T",(si,6',)r=i), H = H{si,...,Sn) &Tldu = Y,ds,^dei. (4.26) 

i 

Definition 12 Let { he affine coordinates in a {n + 1)- dimensional affine space 

A^^^ . A hypersurface in A""*"^ is called a hypersurface of translation if it admits a parame- 
trization of the form 

{y\ . . . - x^{y\ . . . ,2/") = ai(2/i) + . . . + (j = 1, . . . , n + 1). (4.27) 

Theorem 13 ([0) The completely integrable Hamiltonian system ( \i.2di ) admits a second 
Poisson structure, giving its complete integrability if and only if the graph of the Hamiltonian 
function is a hypersurface of translation relative to the affine structure determined by the 
action variables. 
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We present the "only if part of the proof because in what follows the new Poisson 
structure will be needed. 
Proof: 

Assume the {M'^^,uj,H) is a completely integrable system and that the graph of if is a 
hypersuface of translation relative to the action variables (s'), so it has a parametrization 
of the form ( [4.27| ) with = s\ i = 1, . . . ,n and x"'^'^ = H. We can choose the parameters 
(y*) so that the Hamiltonian takes the simple form 

H{y\...,y^)=y' + ... + y\ 

If {tp^ . . . , are the coordinates conjugated to {y^, . . . , y"), we define a second Poisson 
structure by the formula 

M = y2y'^A^. 

^ ay* o^p^ 

One checks easily that the two Poisson structures are compatible, and that the recursion 
operator is given by 

"9 d 
R = Y./i^,®dy^ + ^.®d^')- 

It is now clear from the expression of the Hamiltonian function in the ^/-coordinates that 
Cx„R = 0, so the vector field Xh is bi-Hamiltonian. 



Thus a completely Hamiltonian system whose Hamiltonian's graph is a hypersurface 
of translation relative to the affine structure determined by the action variables, is bi- 
Hamiltonian 

Xh = kl{dH) = K\{dH^) 

n 

withi:fo = EM2/')- 

i=l 

The (1, l)-tensor 

" (9 d 



satisfies 



jCxhP = Id and CpxR = 

n o 



so we can apply Theorem ^ and conclude that the vector fields = R FX = ^ y^cp 

i=l 

n 

and the function G = '^y'^Lp^ allows us to define the sequence of functions Gj = YiiG) {i 

i=l 
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0, 1, . . .), such that 

X,{G)=tHi, 
{X,{G,),Gk} = {t+j)X,{Gk+,) 

and 

{Gi, Gj}o = (i — j)Gi^j. 

Thus, for each z G N, the integrals of motion of Xi, L\ - = Xi{Gk))Gj — Xi{Gj)Gk and 
Xi{Gj), close quadratically under the Poisson bracket defined by Aq. 

Example 14 (The Rational Calogero-Moser System) 

The rational Calogero-Moser system is an integrable Hamiltonian system defined by 

^ = E(y + |-E(* (4-28) 

It admits the pair of matrices (L, M), 

Lij = P^Si + gV^iq^ - q,)-\l - 6{), (4.29) 



M,, = ^?v^E(* - - aV^iq. - qjr'il - Si) (4.30) 

as a Lax pair. 

The Hamiltonian vector field with respect to the canonical Poisson structure in R^"^ is 

This system is completely integrable when we consider the sequence of integrals of motion 
Fi = Tr{L'),i = l,...n. 



Moreover, following [jl2|, if we consider the functions = Tr(QU ), which provide the 



algebraic linearization of the system , |^ , the Hamiltonian vector field becomes 

(4.32) 
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and we may define the following compatible Poisson tensors 

The system u = Xi{u) is bi-Hamiltonian with respect to these Poisson structures and 
the bi-Hamiltonian sequence of integrals of motion is 

Alidh,) = A{{dh,,i), J =0,1,... (4.35) 

where 

= ln(Fi...F„), 

= ^(f^f'-^^^^^o-y^' = T.iF^y^'^ J = 0, 1, ... . (4.36) 

Notice that Xi = Ao(o?/ii) = A\{dho) and if we define the sequence of Hamiltonian vector 
fields 

X, = (A»A«-^)-iX, = JZi^.TT^ ^ = 1, 2, . . . , (4.37) 



the following relation holds 

X,(£GkFu) = {i + l)K (4.38) 

k 

Now consider the (1, l)-tensor 
that satisfies 

CpxR = R and CxP = Id 
and define the sequence of master symmetries 

= (A»A«~^)*PXi. (4.39) 

Considering the functions gi = Zi{J2k{GkFk)) = J2k Fl^^Gk we have 

k 

{Xi{gj),gk}o = {i+ 3 + l)^i(5'i+fc) 
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and 

{9i,9j}o = {i- i)9i+3- 

Now the Propositions ^ and ^ ensure that for each Xj, z G W, the integrals Xi{gj) and 
L\ j = Xi{gk)gj — Xi{gj)gk close quadratically under {., .}o. 



Example 15 (The Goldfish System) The Goldfish system was first introduced by Calo- 
gero [m and is the Hamiltonian system in (IR^", {qi,Pi)) defined hj X = Aq^cIH), where Aq 
is the canonical Poisson tensor and 

H = T^ e"^" (A AO) 

with gi arbitrary smooth functions, each one depending only on the corresponding coordinate 
Qi and a an arbitrary constant. 



Defining the Nijenhuis tensor 



and the function 



R = qi{-rr- ®dqi + — ® dpi) (4.41) 

OQi OPi 



^ = E n"Tf''ll ^°" ^^ = n'^. (4-42) 



this system becomes completely integrable and quasi-bihamiltonian [|T^]. The functions q, 
coefficients of the minimal polynomial of R (g" + Y.iCiq^~^ = Y\!i=i{(l — Qi)), together with 
the integrals of motion 

linearize algebraically the system because = aF^. 

In the coordinates {ci,Fi) the system has a simple bi-Hamiltonian structure defined by 
the compatible Poisson tensors 

and the Hamiltonians 

Ho = F, + ... + F^, H, = ]^{Fl + ... + Fl). (4.45) 
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Similarly to last example, the (1, l)-tensor 

n 

and the function G = ^ FjCj allow us to define the integrals of motion that close quadrati- 

1=1 

cally. 
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